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Abstract. We study the surface brightness profiles of a sample of 25 distant (0.3 < z < 0.83) hot (kT > 3.5 keV) 
clusters, observed with ROSAT, with published temperatures from ASCA. For both open and flat cosmological 
models, the derived emission measure profiles are scaled according to the self-similar model of cluster formation. 
We use the standard scaling relations of cluster properties w ith r edshift and temperature, with the empirical 



slope of the M gas -T relation derived by Neumann & Arnaud (2001). Using a \ 2 test, we perform a quantitative 
comparison of the scaled emission measure profiles of distant clusters with a local reference profile derived from 
the sample of 15 hot nearby clusters compiled by Neumann & Arnaud (1999), which were found to obey self- 
similarity. This comparison allows us to both check the validity of the self-similar model across the redshift range 
0.04 — 0.8, and to constrain the cosmological parameters. 

For a low density flat universe, the scaled distant cluster data were found to be consistent, both in shape and 
normalisation, with the local reference profile. It indicates that hot clusters constitute a homologous family up to 
high redshifts, and gives support to the standard picture of structure formation for the dark matter component. 
Because of the intrinsic regularity in the hot cluster population, the scaled profiles can be used as distance 
indicators, the correct cosmology being the one for which the various profiles at different redshifts coincide. 
The intrinsic limitations of the method, in particular possible systematic errors and biases related to the model 
uncertainties, are discussed. Using the standard evolution model, the present data allow us to put a tight constraint 
on fio for a flat Universe: fio = 0.40lo'i2 a ^ ^0% confidence level (statistical errors only). The critical model 
(Qo = 1) was excluded at the 98% confidence level. Consistently, the observed evolution of the normalisation 
of the Lx~T relation was found to comply with the self-similar model for Qq — 0.4, A = 0.6. The constraint 
derived on Qq is in remarkable agreement with the constraint obtained from luminosity distances to SNI or from 
combined analysis of the power spectrum of the 2dF galaxy redshift Survey and the Cosmic Microwave Background 
anisotropies. 
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1. Introduction 

In the simplest models of structure formation, purely 
based on gravitation, galaxy clusters constitute a ho- 
mologous family. Clusters are self-similar in shape, and 
predictable scaling laws relate each physical property to 
the cluster total mass M and redshift z (Kaiser 



1986 



Navarro et al. |1997 |; Teyssier et al. 1997|; Eke et al. 
Bryan & Norman 



1998). Self-similarity applies to both 



the dark matter component and the hot X-ray emitting 
intra-cluster medium (ICM). 
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From the observation of the ICM, we do see regu- 
larity in the local (z < 0.1) population of clusters, like 
strong correlations between luminosity, gas mass, total 
mass, size and temperature T (Mohr et al. 1997; Allen & 
Fabian |1998| ; Marke vitch [199$ A rnaud & Evrard [1999 ; 



Horner et al. 1999; Mohr et al. 1999 Vikhlinin, et al 



1999; Ncvalainen et al. 2000; Finogucnov et al. 2001; Xu 



1998 ; et al. 2001). Furthermore, there is strong indication of 



a universal shape for the density and temperature pro- 
files of hot (kT > 4 keV) clusters, beyond the cooling 
flow region (Markevitch et al. 1998; Neumann & Arnaud 



1999, 2001; Vikhlinin, et al. 1999; Irwin & Bregman 2000 



Arnaud 2001). However, clusters also deviate from the 



mar- simplest self-similar model. The most remarkable devia- 
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tion is the slope of the luminosity-temperature (Lx~T ) 
relation, which is steeper than predicted. In a recent study, 
Neumann & Arnaud (2001) showed that a steepening of 



the M gas -T relation (Af gas oc T 1 - 94 , instead of the stan- 
dard relation M gas oc T 15 ) can explain the observed Lx~T 
relation in the hot temperature domain, and account for 
the scaling properties of the normalisation of the emission 
measure profiles of hot clusters. Similar steepening was 
derived from direct studies of the M„ as ~T relation inde- 



pendently carried out (Mohr et al. 1999; Vikhlinin, et al 



199E) and is also consistent with the observed slope of the 



isoph otal size temperature (S\-T ) relation (Mohr et al. 



19971) 



Several physical processes have been suggested to ex- 
plain the departure from the simplest self-similar model. 
Pre-heating by early galactic winds, has been proposed to 
explain the steepening of the Lx~T relation (e.g. Kaiser 



1991 ; Evrard & Henry 1991 ), alth ough other effec ts lik e 
AGN heating (e.g. Valageas & Silk |l998fc Wu et al. |2000D , 
radiative cooling (Pearce et al. 200C ; Muawong et al. 



2001) or variation of the galaxy formation efficiency with 



system mass (Bryan 1998) might also play a role. Further 
evidence of the importance of non-gravitational processes 
is provided by the excess of entropy (the "entropy floor" ) 
in poor clusters (Ponman et al. 1999| ; Lloyd-Davis et al. 
2000| ). Recent numerical simulations (Bialek et al. 09991) 
including pre-heating, with an initial entropy level con- 
sistent with this observed entropy floor, do predict a 
steepening of the Lx~T , M gas -T and S\-T relations, 
consistent with the observations quoted above (see also 
Loewenstein 2000; Tozzi & Norman 2001; Brighenti & 
Mathews 2001; Borgani et al. 2001). However, it is unclear 
if such a scenario is also consistent with the level of self- 
similarity in shape observed in hot clusters. Although it is 
predicted that cool clusters should have a more ex tende d 
atmosphere than hot clusters (e.g. Tozzi & Norman 2001 ), 
to our knowledge no detailed study on the relationship be- 
tween internal shape and cluster temperature, specifically 
for relatively hot clusters, has been carried out so far. 

The evolution of cluster X-ray properties is an essen- 
tial piece of information to reconstruct the physics of the 
formation processes for the gas component and can also 
be used as a cosmological test. Models with pre-heating 
predict an absence of evolution in the Lx~T and M gas -T 
relations, at least up to z ~ 0.5 (e.g. Bialek et al. 1999| ). 
There is some indication, based on a few massive clus- 
ters, that the Lx~T relation is evolving weakly, if at all 
(Mus hotzky & Scharf |l997j; Sad at et al. |1998|; D onahue 
Reichart et al. |1999| ; Schindle r^l999| ; Fairley 



199E 



et al. 

et al. 200C). Several groups (Sadat et al. 1998; Reichart 
Fairley et al. 2000D quantified the evolution 
of the normalisation of the Lx~T relation, assuming it 
varies as (1 + z) v . For a critical density Universe, they 
found rj values significantly smaller than the theoretical 
prediction r\ = 1.5 in the self-similar model and consistent 
with no-evolution. However, the luminosity estimates de- 
pend on the assumed cosmological parameters and so does 
the constraint on the evolution parameter (Fairley et al. 



|200C ; Reichart et al. 1999). The evolution of other scaling 
laws like the gas or total mass temperature relation are 
even more poorly known (Schindler 1999; Matsumoto et 
al. pOOOj). 



Using non-evolving physical properties of clusters as 
distance indicators can provide interesting constraints on 
cosmological parameters, such as the density parameter, 
f2oi and the cosmological constant, A. In this context, 



the gas mass fraction has been considered by Pen (1997), 
although present constrains are poor (Rines et al. |1999; 



Ettori & Fabian |1999[ ). Recently, Mohr et al. Q200C| ) mea- 
sured the S\ T relation for a sample of intermediate red- 
shift clusters, 0.2 < z < 0.55. Using standard cluster evo- 
lution models, they argue that this relation should not 
evolve with redshift. They did find that the intermediate 
redshift data are consistent with the local relation and 
were able to rule out a critical density Universe. 

With the present study, we aim at a better understand- 
ing of the evolution of the scaling and structural properties 
of hot clusters with redshift. Furthermore, we show that 
strong constraints on the cosmological parameters can be 
drawn, based on the cluster scaling properties. 

We perform for the first time a systematic study of 
the X-ray surface brightness profiles of distant (0.3 < 
z < 0.83) hot (kT > 3.5keV) clusters, measured with the 
ROSAT satellite. This sample is combined to the sam- 
ple of loc al (z ~ 0.05) clusters, presented in Neumann & 
Arnaud (1999). The surface brightness profile is directly 
related to the emission measure profile (or equivalently to 
the gas density profile). Comparing the profiles of clusters 
at different redshifts and temperatures obviously provides 
more information than simply considering global quanti- 
ties such as the total X-ray luminosity or punctual quan- 
tities like the isophotal radius. With the present study, 
we wish to address the following issues i) Do hot clus- 
ters remain self-similar in shape up to high redshift? ii) 
How do the scaling properties of the profiles with redshift 
compare quantitatively with the theoretical expectations 
of the self-similar model? iii) What constraints can we put 
on the cosmological parameters from these data? iv) Is the 
evolution of the Lx~T relation really inconsistent with a 
self-similar model? 

The paper is organized as follows. In Section 2, we 
present the cluster sample and the data analysis per- 
formed to derive the surface brightness profiles and then 
the emission measure profiles. In Section 3 we derive how 
the emission measure profiles should scale with redshift, 
depending on cosmological parameters, for the self-similar 
model of cluster formation. In Section 4, we derive the cor- 
responding scaled emission measure profiles for our clus- 
ter sample, that we use, in Section 5, to test the self- 
similar model and constrain the cosmological parameters. 
In Section 6, we study the the Lx~T relation. In Section 
7 we discuss our results and Section 8 contains our con- 
clusions. 

The present time Hubble constant in units of 50 
km/s/Mpc is noted I150 in the following. The data analysis 
is done with hso = 1. 
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Table 1. Basic data for the 25 distant clusters in the sample. 
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^det 
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0.308 
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12 -l-0.5 


5.52 
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(P) 


25.2 ± 0.5 
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4.88 


7.1 ± 0.1 


CLG 0016+1609 


0.541 
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7.29 


4 


43 


(P) 
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5.21 
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0.373 


n ~+0.6 
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1.75 


7 


32 
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2.5 ± 0.1 


2.0 


2.15 


2.3 ± 0.1 


CL 0302.7+1658 


0.424 




1.08 


4 


34 


(h) 


0.72 ± 0.07 


1.0 


1.64 


1.0 ± 0.1 


MS 0353.6-3642 


0.32 


n r + l 
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1.43 


4 


22 
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2.7 


2.45 
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0.55 


1 „ O+0.9 
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4 


16 
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4.5 


4.98 


6.1 ± 0.2 


"TV/TO ao 1 -i r i i /"i o i 

MS 0811.6+6301 


0.312 


4 - 9 -o.6 


0.570 


4 
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(h) 


0.96 ± 0.05 


1.3 


2.02 


0.58 ± 0.03 


MS 1008.1-1224 


0.301 


o n+1 2 


1.84 


4 


69 


(h) 


2.6 ± 0.1 


2.3 


1.72 


2.0 ± 0.1 


A 959 


0.353 


r*+l.l 

7-OIo.g 


2.30 


6 


16 


(P) 


7.2 ± 0.2 


5.5 


5.40 


1.81 ± 0.06 
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0.83 


„ _ _i_2 l 


4.39 


6 
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2.30 
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1 1 ,+0.8 
n - 4 -0.6 
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8 


8.6 
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5.5 


4.56 
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MS 1137.5+6625 
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c 7+1.3 


1.62 


2 


99 


(h) 
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1.0 


2.80 


2.0 ± 0.1 


MS 1224.7+2007 


0.327 


. - +0.7 
ffc " L — 0.6 


0.690 


4 


41 


(h) 


0.82 ± 0.07 


1.0 


2.19 


0.58 ± 0.05 


Mb 1241.5+1710 


0.54 




2.26 


4 


31 
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1.3 ± 0.1 


1.2 


2.48 
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A 1 TOO 
1\ L ( ZZ 


U.oZlD 


r Q+0.3 
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(h) 


9 10 






1.0 ± U.l 


RX J1347.5-1145 


0.451 


O Q+0-7 
a '°-0.6 


21.0 


9 


36 


00 


12.7 + 0.3 


3.0 


1.96 


17.4 + 0.4 


Zwcl 1358+6245 


0.328 


6 9 +0 - 5 

°- y -0.5 


2.14 


4 


23 


(P) 


8.0 + 0.2 


2.7 


5.20 


1.95 + 0.05 


3C295 


0.46 




1.90 


6 


29 


(h) 


2.1 + 0.1 


1.0 


2.52 


2.9 + 0.1 


MS 1426.4+0158 


0.32 


6.4±i;° 


0.970 


4 


37 


(h) 


1.4 + 0.1 


1.3 


2.16 


1.09 + 0.08 
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0.318 


10.7+lf 


2.82 


6 


38 


(h) 


4.3 + 0.2 
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2.21 


3.7 + 0.1 


MS 1512.4+3647 


0.372 


o 4 +0.4 
"■>•*— 0.4 
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4 


35 


(h) 


1.8 + 0.1 


2.0 


2.52 


0.86 + 0.07 


MS 1621.5+2640 


0.426 


6 6 +0 - 9 


1.58 


4 


44 


(h) 


1.06 + 0.09 


1.3 


2.13 


1.7 + 0.1 


RX J1716.4+6708 


0.813 


c 7 +1.4 


1.15 


3 
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(h) 


0.37 + 0.03 


1.0 


2.34 


1.5 + 0.1 


MS 2137.3-2353 


0.313 


4 q+0-3 


3.35 


4 


10 


(P) 


16.2 + 0.4 


3.5 


5.06 


2.91 + 0.07 


ACCG 114 
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q 8 +0 {i 


3.25 


1 


23 
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6.5 + 0.3 


3.8 


2.26 
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Notes: The values of the bolometric luminosities 
68% confidence level. 
a The redshifts, z, are taken from NED. 

b 



Lbd, are for Qq — 1 (qo = 0.5) and Ho = 50kms/s/Mpc. All errors are at the 



al 



Referen ces for the tem perat ure and A SCA l uminosities listed c olum n (3) and (4): 
3. Gioia et al. 



1999 



Pierre et al. 



1999 



4. Henry 
19991; 9. Schindler. et al 



2000 



1997 



2001; 6. Mushotzky 



Allen & Fabian 
i Scharf 



1998| ; 2. Do nahue et 
19971; 7. Ota et al. h998|; 8. 



5. Jeltema et al 

The temperature errors published at the 90% confidence level were divided by 
1.65 to estimate the 68% confidence level errors. Luminosities published in the [2 — 10]keV energy band (references 4. and 7.) 
were converted to bolometric luminosity using a MEKAL model with the temperature given column (3). When necessary, the 
published luminosities were corrected for Ho = 50kms/s/Mpc and qo = 0.5 (references 2., 3. and 6.). 
c The letter in parenthesis stands for the ROSAT detector used, (h) for HRI, (p) for PSPC. 



2. The data 

2.1. The cluster sample 

We considered all distant (z > 0.3) clusters observed by 
ROSAT, with published ASCA temperatures. We believe 
our original list was complete with respect to ROSAT pub- 
lic archival data and publications, available at the end of 
1999. We excluded three clusters with no obvious X-ray 



center: the double cluster A851 (Schindler et al. 1998| ), 
the clumpy cluster CI 0500-24 (Schindler & Wambsganss 



1997) and MS 1147.3+1103 (the HRI image shows a very 



flat elliptical morphology in the core, with some evidence 
of bimodality). The derivation of a surface brightness pro- 
file for those clusters would have been arbitrary. We also 
exclu ded CI 2244-0221 and MG 2053.7-0449 (Hattori et al. 



1997) due to the too poor statistical quality of the HRI 



The list of the 25 distant clusters selected is shown 
in Tab.|l], as well as the exposure times and the ROSAT 
detector used. The sample covers a redshift range of 
z = 0.3 — 0.83. We also give in the table the tempera- 
tures and bolometric luminosities, measured with ASCA. 
The only exception is MS1054, for which we list the recent 
Chandra temperature estimate of the main cluster corn- 



data. 



ponent (Jeltema et al. 2001), the western subcluster (see 
Neumann & Arnaud 2000} ) being excluded in our spatial 
analysis below. When several temperature estimates for 
a given cluster were published, we have chosen the most 
recent analysis using the latest ASCA calibrations. The 
various published values were usually consistent. 

To study cluster evolution, we combined this new dis- 
tant cluster sample with the sample considered in our pre- 
vious study of the surface brig htness pro files of nearby 
clusters (Neumann & Arnaud 1999 , 2001 ). This nearby 
cluster sample comprises 15 Abell clusters in the redshift 



4 



Arnaud et al.: The X-ray surface brightness profiles of hot galaxy clusters up to z ~ 0.8 
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Radius (arcmin) Radius (Mpc) 

Fig. 1. The ROSAT X-ray surface brightness (Sx) profiles (left panel) and emission measure (EM) along the line of 
sight profiles (right panel) of all (nearby and distant) clusters in the sample. A straight line is drawn between the data 
points for each cluster for a better visualization of the profiles. The curves are color-coded by redshift. The Sx profiles 
are background subtracted and corrected for vignetting effects. EM is derived from the surface brightness taking into 
account the (1 + z)~ 4 cosmological dimming and the emissivity in the ROSAT band (see Eq. The angular radius 
is converted to physical radius, assuming fl — 1 and h 50 = 1. 



range 0.04 < z < 0.06, which were observed in pointing 
mode with the ROSAT PSPC with a high signal to noise 
ratio and for which accurate temperature meas ureme nts 
exist from the literature (see Neumann & Arnaud 1999 for 
details). 

We emphasize that the study presented here focuses 
on relatively hot clusters, the minimum temperature for 
the nearby and distant cluster samples being 3.7 and 3.4 
keV respectively. 

2.2. Surface brightness profiles 

The surface brightness profile of each cluster, S(0), was 
constructed using the standard procedures described in 
Neumann & Arnaud ( 1999[ ). We only considered photons 
in the energy band 0.5-2.0 keV for the PSPC data and only 
took into account channels 2-10 for the HRI data, in order 
to optimize the signal-to-noise (S/N) ratio. We binned the 
photons into concentric annuli centered on the maximum 
of the X-ray emission with a width of 15" and 10" per 
annulus for the PSPC and HRI data respectively. We cut 
out serendipitous sources in the field of view or cluster 
substructures, if they show up as a local maximum. The 



HRI particle background was subtracted from the HRI 
profiles using the background map construct ed for each 
observation with the method of Snowden et al. fll99Sp . The 
vignetting correction was performed using the ex posure 
maps computed with EXSAS (Zimmermann et al. 1994) 
for the P SPC d ata and with the software developed by 
Snowden ( 1998 ) for the HRI data. The X-ray background 
for each pointing was estimated using vignetting corrected 
data in the outer part of the field of view and subtracted 
from the profile. A 5%(10%) systematic error was added 
quadratically to the statistical error on the PSPC(HRI) 
background level. 

For 5 clusters both HRI and PSPC data were available. 
We found an excellent agreement between the HRI and 
PSPC profiles, except for the most inner radial bin, where 
the effect of the wider PSPC/PSF can be observed. This 
blurring is clearly negligible at larger radii. If available, we 
thus always choose the PSPC data, due to its higher in- 
trinsic sensitivity and lower background level, which allow 
to trace the cluster emission further out. 

To avoid too noisy profiles, we rebinned the data, for 
both the nearby and distant cluster samples, so that the 
variations of S(9) from bin to bin are significantly larger 
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than the corresponding statistical error and thus repre- 
sentative of the cluster shape. Starting form the central 
annulus, we regrouped the data in adjacent annuli so that 
i) at least a S/N ratio of 3cr is reached after background 
subtraction and ii) the width of the annulus, A(#), at ra- 
dius 6, has a size at least 0.15 6. This logarithmic bin- 
ning insures a roughly constant S/N ratio for each bin in 
the outer part of the profiles, where the background can 
still be neglected (the S/N ratio would be constant for a 
/3-model with (3 — 2/3 and no background). The adopted 
rebinning was found to be a good compromise between the 
desired accuracy and a reasonable sampling of the profiles. 

The resulting surface brightness profiles are shown in 
Fig. [l] (left panel) for the distant and nearby clusters. For 
each cluster the data points are connected by a straight 
line to guide the eye. The profiles are plotted up to the 
adopted detection limit of 3cx above background. The cor- 
responding detection radius, as well as the total ROSAT 
count rate, Ci?dot, within this radius is given in Tab. [I] 
for each cluster. In Fig. [l] and Fig. ||, we color-coded the 
profiles according to cluster redshift: blue for nearby clus- 
ters (0.04 < z < 0.06), green for moderately distant clus- 
ters (0.3 < z < 0.40) and red for very distant clusters 
(0.40 < z < 0.83). Note that this redshift subsampling 
of the distant cluster sample is for display only and, un- 
less explicitly stated, is not used in the statistical analysis 
below. 

2.3. Emission measure profiles 

The emission measure along the line of sight at radius 
r, EM(r) = j n c nn dl, can be deduced from the X-ray 
surface brightness, S(8): 



EM(r) 



4 7T (l + z) 4 S(0) 
e(T,z) 

4.81 HT 7 (l + z) 4 ( 

e(T,z) 



S(0) 
ct/s/arcmin 2 
l 



K ct/s/10 10 cm- 5 



cm 



Mpc 



(1) 
(2) 



where d\(z) is the angular distance at redshift z. 

e(T, z) is the emissivity in the considered ROSAT 
band, \E\ — E^], taking into account the interstellar ab- 
sorption and the instrumental spectral response: 



e(T,z) 



S(E)e 



~a(E)N H 



f T ((l + z)E){l + z) 2 dE (3) 



El 



where S(E) is the detector effective area at energy E, <j(E) 
the absorption cross section, Nn the hydrogen column 
density along the line of sight, fx((l + z)E) the emissivity 
in photons cm 3 /s/keV at energy (1 + z)E for a plasma of 
temperature T. It was computed for each clus ter using a 
redshiftc d the rmal emission model (Mewe et al. 1985 , 1986 ; 
Kaastra [l99^ ; Lieda hl et a l. |l995| ), the ROSAT response 
(Zimmermann et al. 1994) and the Nn value estimated 



with the w3nh tools available at HEARSAC (Dickey & 
Lockman |199C ). The emissivity, e(T, z), depends weakly 
on cluster temperature and redshift in the ROSAT band 
(Tab. §. 

The derived EM(r) profiles are shown in the right 
panel of Fig. [l]for a critical density Universe (f2o = 1, A = 
0.). As will be discussed later, one can already note that 
the distant clusters appear brighter than the nearby clus- 
ters. 



3. Theoretical scaling laws 

In our derivation of the theoretical emission measure pro- 
files, as a function of redshift and cosmological parameters, 
we will consider both a flat Universe (Slo + A = 1) and an 
open Universe (flo < 1,A = 0). The matter density pa- 
rameter at redshift z is noted f2 z ; £! z = £!o(l + z) 3 /E 2 (z), 
where E 2 {z) = Q (l + z) 3 + (1 - O - A)(l + z) 2 + A. 



3.1. The self-similar model 

The simplest self-similar model (e.g. Bryan & Norman 



1998; Eke et al. 1998) assumes that i) at a given redshift 



the relaxed virialized portion of clusters corresponds to a 
fixed density contrast as compared to the critical density 
of the Universe at that redshift ii) the internal structure 
of clusters of different mass and z are similar. 

The virial mass My and radius Ry then scale with 
redshift and temperature via the well known relations: 



My = 2.835 10 15 /4 /2 AJ 1 / 2 (1 + z) 



kT x 3/2 



-3/2 



h^ Mr. 



Ry — 



with 
A z 



10 keV 
3.80 (3 X 4 2 A- 1 / 2 (1 + z)- 3 / 



50 m & 
-3/2 



kT 



10 keV 



1/2 



(A c (ft z ,A)f} )/(18A) 



(4) 

(5) 
(G) 



where A c (f2 z , A) is the density contrast (a function of f?o 
and A) and /3t is the normalisation of the virial relation, 
GMy/2Ry = /3 T kT. 

The My-T and Ry-T relations depend on the cos- 
mological parameters through the factor A c (f2 z , A)f2 /^z- 
This factor is constant with redshift and equal to 187r 2 
for a critical density Universe. Analytical approximations 
of A c (f2 z ,A), derived from the top-hat spherical collapse 
model assuming that clusters have just virialized, are given 
in Bryan & Norman ( |1998| ): 



A c (O z ,A) 
A c (O z ,A) 
with w 



= 187T 



18 
18 

a 



= 18vr 



- 60w 
-82w 



1 



- 32w 2 for Q < 1,A = 
■ 39w 2 for n Q + A = 1 



(7) 



As we consider lower and lower values of the density pa- 
rameter SIq, the assumption of recent cluster formation 
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is less and less valid and in principle, the difference be- 
tween the observing time and the time of collapse has to 
be taken into account (e.g. Voit & Donahue 1998). If the 
effective formation epoch of a cluster of a given mass is 
earlier than the observing time, when the Universe was 
denser, the actual cluster temperature is underestimated 
by the recent formation approximation. This effect is ex- 
pected to increase with decreasing Qq and mass of the 
system. Estimating accurately the impact on the My-T 
relation (mean relation and scatter) is however not triv- 
ial, because it requires a precise modelling of cluster for- 
mation history, the growth of clusters by continuous ac- 
cretion and m erger events, and the comp lex physics of 
the ICM (Voit |2000l Afshordi & Cen |200 j) . However, re- 
cent analysis of the My-T relation derived from numerical 
simulations suggest that the effect of formation redshift 
is negligible, at least when considering measured X-ray 
temperatures (Mathiesen 1998 ). We will thus neglect this 
effect here and use the above equations estimated at the 
observed cluster redshift. 

The constant /3t depends on the cluster internal struc- 
ture. Its value can be determined from numerical simula- 
tions. The various results agree within typically ±20%, 
with no obvious dependence on cosmological parameters 
(Henry 2000). As in our previous work (Neumann & 
Arnaud 1999;2001), we will adopt the normalisation of 
Evrard et al. (1996), p\ — 1.05. Note that our results do 
not depend on the exact value of p\. 



Assuming a standard /3-modcl with (3 = 2/3 and a scaled 
core radius of x c = 0.123, which fits well the scaled profiles 
of nearby clusters (Neumann & Arnaud 1999), gives Q n = 
69.4. 

The scaling law for the central emission measure can 
now be derived from Eq. |8|JTc| , |TT| and |l^ and the My-T 
and Ry-T relations (Eq. |4,|5j), assuming that all clusters 
have the same gas mass fraction / gas : 

1/2 / £ \ 2 



EMo = 4.1 10" 



1.05 



x A 3 / 2 (1 + zf' 2 



/gas \ f Qn 

YT \ 1/2 
10 keV J 



h 3 



5 Mf£<3) 



This assumes that the gas mass scales as the total mass, i.e 
M gas ex T . The corresponding emission measure profiles 
can thus be written: 



EM(r) = EM hm{r/Ry) 
where /em(#) is the dimensionless function: 
/»(«) 



/em (x) 



y^ 2 



=d(u 2 



(14) 



(15) 



As can be seen from Eq. |^ and Eq. 13, the virial radius 
decreases with redshift while the central emission measure 
increases. Clusters of a given mass are denser at high red- 
shift, following the evolution of the Universe mean density. 
We thus expect that clusters of given temperature appear 
smaller and brighter with increasing redshift. 



3.2. The theoretical emission measure profiles 

The central emission measure along the line of sight is 
related to the electron density profile of the gas, n c (r), 
via: 



rRv 

EM = 2 (n H /n e ) / n c (r) 2 dr 
Jo 



whereas the gas mass is given by: 



Mgas = A*' m p («h/«-c) 47T 



n c (r) r 2 dr 



(8) 
(9) 

(10) 



where m p is the proton mass, // = 1.347 and nn/fi c = 
0.852, for an ionized plasma with a metallicity of 0.3 solar 
value. In self-similar models, which we consider here, the 
density profile can be written: 



n c (r) = n c (0) f n (x) ;x = r/Ry 



(11) 



where x is the radius scaled to the virial radius and / n is a 
universal function, the same for all clusters. By combining 
the above equations, EM varies as EMo oc Q n M 2 /R 5 



v 



where we have introduced a constant form factor Q n , 
which only depends on the cluster's 'universal' shape: 



Qn = 



x 2 dx 



(12) 



3.3. The Lx~T relation and the empirical scaling law 



The bolometric cluster luminosity is given by: 
L x = A(T) / EM(r)2Trrdr 



(16) 



where the cooling function, A(T) varies as A(T) cx T 1 / 2 . 
For the standard self-similar model described above, the 
bolometric luminosity follows the well known scaling rela- 
tion: 



L x oc A, 1 / 2 (1 + zf' 2 T 2 



(17) 



which is inconsistent with the slope of the observed local 
L x -T relation, a - 2.88 (e.g Arnaud & Evrard |l999| ) . 

As already mentioned in the introduction, we found 
evidence for a steepening of the M gas -T relation for hot 
clusters (Neumann & Arnaud 2001 ), in our previous study 
of the nearby cluster sample considered here. A gas mass 
varying as M gas oc T 194 , instead of M gas oc T 15 , can both 
explain the observed L-^-T relation and significantly re- 
duce the scatter in the scaled emission measure profiles, 
when compared to the standard scaling. In that case, the 
emission measure scales with temperature as EM oc T 138 , 
instead of T 5 . We will also consider this empirical scal- 
ing law in the following section. The dependence of the 
normalisation on redshift and cosmological parameters re- 
mains a priori unchanged and we will assume that the em- 
pirical slope of the relation does not evolve with redshift, 
which is the simplest assumption. 
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Fig. 2. The scaled emission measure along the line of sight as a function of scaled radius for the 40 clusters of the 
sample in two cosmological models. Top panels: Oq = l.,A = 0. Bottom panels: Hq — 0.3, A = 0.7. The data points 
for each cluster are connected by a straight line for a better visualization of the profiles. The curves are color-coded 
by redshift. The radius is scaled by the virial radius, Ry, computed using the theoretical scaling law (Eq. ||) with the 
measured temperature and redshift of the clusters. The emission measure profiles in the left panels have been scaled by 
(A c (fi Z) A)fi /^z) 3/2 (l + z) 9 / 2 T 1 / 2 , according to the standard self-similar model (Eq. Q, with M gas oc T 15 . Note the 
remarkable similarity of the profiles at radii larger than ~ 0.1i?v- The scaling procedure has significantly decreased the 
dispersion between the emission measure profiles (compare with the right panel of Fig|l]). The right panels correspond 
to the e mpiric al scaling law assuming Af gas oc T 1 - 94 , consistent with the slope of the local Ly^-T relation (Neumann & 
Arnaud 2001). The dispersion is decreased even more as compared to the standard scaling. 
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4. Scaled emission measure profiles 

4.1. Scaling procedure 

As in our previous studies, we scaled the emission mea- 
sure profiles so that they would lie on top of one another if 
obeying self-similarity The scaled emission measure pro- 
files, corresponding to the standard scaling with z and T 
given Eq. n~3], are thus definedP] as: 



-1/2 



X 

r 
Ry 



EM(r) 



4.10 10- 6 cm~ 6 Mpc 



(18) 



where Ry is defined by Eq. [| and the emission measure is 
derived from the surface brightness via Eq. |l|. 

To introduce the empirical EM-T scaling relation 
(EM cx T 1 ' 38 ), we simply have to introduce a corrective 
multiplicative factor of cx J 1-0 - 88 to the previous equation: 



EM(x) = A" 3 / 2 (l + z)- 9 / 2 
EM(r) 



kT 



10 kcV 



-1.38 



6.0 10" 



>Mpc 



(19) 



For convenience, the corrective factor has been arbitrarily 
normalized to 1 for a temperature equal to 6.5 keV, which 
is the mean temperature of the sample. 

The scaled profiles depend on the assumed cosmolog- 
ical parameters, via the angular distance d\(z) used to 
convert angular radii to physical radii and via the factor 
A z = (A c (f2 z , A)i7o)/(187r 2 f2 z ) appearing in the normali- 
sation of the profiles and of the Ry-T relation. The varia- 
tion with redshift of both quantities depends on fio and A. 
Therefore, if the self-similar evolution model is valid, the 
scaled profiles of clusters observed at various redshifts will 
coincide, but only for the correct cosmological parameters. 

The scaled profiles can thus be used both to check the 
validity of the self-similar model and to put constraints on 
the cosmological parameters, J7o and A. This is described 
in detail in Sect. || and further discussed in Sect. 0. To do 
so, we will use some general properties of the profiles that 
we outline below. 

4.2. Scaled Profiles 

The scaled profiles are shown in Fig. || for two cosmologi- 
cal models, a critical density Universe (f2o — 1) and a flat 
model with Slo = 0.3 and A = 0.7. At first sight, distant 
clusters appear remarkably similar to nearby clusters, once 
the profiles are scaled. As expected, the difference between 



1 Note that the normalisation has been set so that the central 
value would be unity for a /3-model with — 2/3 and x c = 
0.123, a 20% gas mass fraction (see Eq. y_3|) and hso = 1. All 
these factors are common to all profiles and their exact value 
does not matter to check self-similarity 



the scaled profiles, however, depends on the assumed cos- 
mological parameters. The profiles plotted in the top and 
bottom panels of Fig. || are clearly different, in particu- 
lar in the relative position of the clusters for the different 
redshift ranges. 

One further notes the large scatter in the scaled profiles 
in the cluster core (r < 0.1i?v) and the remarkably com- 
mon shape above typically 0.1 — 0.2 Ry. This was already 
noted for nearby clusters (Neumann & Arnaud 1999) and 



clearly also holds for distant clusters. The large scatter in 
the core is likely to be due to cooling flows of various sizes. 
For instance the four clusters with the highest central 
value, MS 2137.3-2353, RXJ 1347.5-1145, MS 1512 and 
3C295, are known, or suspected, to host massive cooling 
flows. The mass accretion rate, estimated using standard 
cooling flow models, is as high as 3000Mf^/yrs for RXJ 



1347.5-1145 (Schindler etal. |l997| and 500 - 9OOM /yrs 
for 3C295 (Neumann 1999 ). The presence of a massive 
cooling flow in MS 1512.4+3647 is indicated by the de- 
tect ion o f luminous extended H a emission (Donahue et 
al. 1992). No detailed cooling flow analysis is available 



for MS 2137.3-2353. However, we note that the estimated 
cooling time for this cluste r is about 1/10 of the Hubble 
time (Allen & Fabian 1998 ) and a clear drop of tempera- 
ture is observed in the center, similar to the one observed 
in RXJ 1347.5-1145 (Allen et al. |200l| ). 

A first quantitative check of similarity beyond the core 
can be made by looking at the dispersion among the pro- 
files at a given radius, for the whole cluster sample. The 
surface brightnesses are measured at discrete values of the 
angular radii 0. To compute the mean value and dispersion 
of the profiles at any physical or scaled radius, a contin- 
uous profile was generated for each cluster using a loga- 
rithmic interpolation of the data. 

The scaling procedure always significantly reduces the 
differences among the profiles, as can already be seen by 
comparing Fig. and Fig. |. This is a first indication 
that clusters obey scaling laws up to high redshift. Let us 
for instance consider the standard scaling with Qq = 1. 
The relative bi- weight dispersion of the emission measure 
profiles at a given radius is ~ 100% between 0.5 and 1 
Mpc, whereas the dispersion drops to ~ 40 — 45% for the 
scaled profiles between 0.2 and 0.5Ry. 

Furthermore, in the same range of radii, the scatter is 
further decreased to ~ 35 — 40% for the empirical scal- 
ing relation (EM cx T 138 , right panel). This decrease 
is slightly more pronounced in a low Qq Universe. The 
improvement is not as spectacular as for the local sam- 
ple alone (a factor of 2 decrease of the scatter). However, 
this additional T~°- 88 scaling factor introduces additional 
noise due to the uncertainties on the temperatures. These 
errors are particularly large for the distant cluster sam- 
ple. The fact that there is still an improvement, in spite of 



2 Each data point is actually the mean surface brightness in 
the radial bin considered and not the surface brightness at the 
center of the bin as assumed here and in the following. We 
checked that the difference is negligible. 
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this additional noise, suggests that the empirical EM-T 
scaling relation fits better the cluster properties than the 
standard case, over the redshift range z — 0.04 — 0.8. We 
will thus adopt this empirical scaling relation in the fol- 
lowing. 

5. Test of self-similarity and constraints on the 
cosmological parameters 

5.1. Method 

Our aim is to check, in a quantitative way, the validity 
of the self-similar model, and set constraints on the cos- 
mological parameters. For that purpose, we need a better 
statistical estimator than the calculated dispersion of the 
profiles at a given radius, which we used in the previous 
section. The relative dispersion is not a global estimator 
and furthermore does not take into account measurement 
errors. 

We first derived, for each set of cosmological parame- 
ters, a scaled reference profile, and an estimate of the in- 
trinsic scatter around it, using the nearby cluster sample 
data. To do so, we estimated, at any given scaled radius, 
the mean value of the different scaled EM profiles, to- 
gether with the corresponding standard deviation, at that 
specific radius. We computed this reference profile up to 
the radius for which at least two nearby cluster profiles 
are still available. Note that measurement errors, which 
are much less than for the distant cluster sample, can still 
contribute to the scatter. Analytical fits of the reference 
scaled profiles (for open and flat Universes) are given in 
Appendix |a|. 

We then considered the set of data points for the dis- 
tant clustersample. Each data point is the scaled emission 
measure EM^j of cluster j, measured at the scaled ra- 
dius Xij, with corresponding 1 a errors. The error on the 
temperature contributes to both the error on a^j and on 
EMij, while the error on the surface brightness obviously 
only contributes to the later quantity. These data points 
are compared to the corresponding reference profile in the 
left panel of Fig. [| for a critical density Universe. 

If the self-similar model is valid, the distant cluster 
data points after scaling must be consistent with the ref- 
erence profile, within the errors. We thus computed the x 2 
value of the distant cluster data about the reference curve, 
for each cosmological model. This x 2 value can be used to 
assess in the standard way the validity of the underlying 
self-similar model and to constrain the cosmological pa- 
rameters, considered as free parameters of the model. 

The x 2 computation is not straightforward, because 
there are non negligible errors on both variables x and 
EM, these errors are correlated, and the reference curve 
is not linear. Furthermore, we have to take into account 
the existence of intrinsic scatter. The computation of x 2 
is detailed in Appendix [b|. 

Another technical issue is the choice of data points in- 
cluded in the computation of the overall \ 2 - First, and 
obviously, only points for which there is a corresponding 
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Fig. 3. x 2 value of the distant cluster data about the mean 
scaled emission measure profile of nearby clusters, as a 
function of £Iq. Full line: flat Universe, Oo + A = 1. Dotted 
line: Open Universe, A = 0. 



reference value from the local sample can be included. In 
practice, very few points are excluded this way, since, for 
every cosmological model, the distant clusters are usu- 
ally traced up to smaller scaled radii when compared to 
nearby clusters. Furthermore, as discussed in the previ- 
ous section, the core properties are clearly dominated by 
different physics. It is thus better to exclude the central 
points to check self-similarity at large radii and to con- 
strain in a more significant way the cosmological parame- 
ters. For that purpose, including the central points would 
be equivalent to add extra noise. We thus considered a 
fixed number of points, N p , defined as the N p most dis- 
tant from the center in scaled coordinates. Although the 
relative position of the points depends somewhat on the 
cosmological parameters, essentially the same data set is 
compared to the reference curve in all casesf]. We both 
considered N p = 80 and N p = 150 corresponding respec- 
tively to a minimum scaled radius x = 0.2 and x = 0.1 for 
Qo = l- 

The variation of the x 2 value with Oo is plotted in 
Fig. H for a flat Universe (Oo + A = 1) and an open 
Universe (A = 0) for N p = 80. The individual x values 
for each data point are plotted on the bottom panels of 
Fig. [| for a critical Universe (left panel) and for the flat 

3 This would not be the case, if we had considered a fixed 
region in terms of scaled radii. The absolute position of the 
profiles in the log-log plane is very sensitive on the cosmology 
as can be seen by comparing the left and right panels of Fig. ^| 
A given angular radius corresponds to a smaller scaled radius 
for a smaller value of Qo- This would have introduced bias 
in the x 2 estimate, with more and more points from the core 
included in the sample as SIq increases. 



10 



Arnaud et al.: The X-ray surface brightness profiles of hot galaxy clusters up to z ~ 0.8 




Fig. 4. Comparison between the scaled emission measure profiles of the distant clusters (data points) and the mean 
scaled profile determined from the nearby cluster sample (z = 0.04— 0.06, full line) for two cosmological models. The 
profiles are scaled assuming M gas oc T 1 - 94 . The data for the distant clusters are the same as shown in the right panels 
of Fig. H but each data point is now displayed individually with la error bars on both scaled variables. Open and 
filled points correspond to clusters in the redshift range z = 0.3 — 0.4 and z — 0.40 — 0.83, respectively. The dotted 
lines correspond to the mean profile of nearby clusters, taken as a reference profile, plus or minus the corresponding 
standard deviation. The bottom panels show the \ value of individual points about the reference profile, taking into 
account the correlation between the errors and the dispersion around the reference profile. Left panel: Results for 
Qq = 1, A = 0. The distant cluster data are not consistent with the nearby cluster data for a critical density Universe. 
Most of the points lie significantly below the reference profile. Right panel: Results for fio = 0.4 and A = 0.6. Note the 
good agreement between the scaled emission measure profiles of distant clusters and the mean scaled profile of nearby 
clusters. There is also no systematic variation of the \ value with radius, an additional indication of self-similarity in 
shape. 



0.6, with the lowest x 2 value N p 



model, ilo = 0.4 and A 
(right panel). 

5.2. Results 



For Slo = 1, most of the scaled distant cluster data points 
fall significantly below the reference curve (Fig. [|, left 
panel) with an overall x 2 value of 108 for N p = 80 
points. Our data thus allow to exclude the critical den- 
sity Universe model at the 98% confidence level. 

Allowing f^o to vary, and considering first a flat 
Universe, we find an excellent agreement of the distant 
cluster data with the reference profile for a low density 
Universe (Fig. ||, right panel). It indicates that hot clusters 
constitute a homologous family up to high redshifts and 
strongly supports the underlying self-similar model. The 
smallest x 2 1S reached for Q = 0.4 (A = 0.6), with x 2 = 77 
(reduced x 2 ~ !)■ Furthermore, no systematic variation of 
the x value with radius is observed, an additional indica- 
tion of the self-similarity in shape of the profiles (bottom 
right panel of Fig ||) . Interestingly a strong constraint can 
be set on fio: Slo = 0.4+q'i2, with errors given at the 90% 
confidence level (corresponding to a A% 2 = 2.7). This re- 
sult is not sensitive to the number of points considered. For 



150, we obtain similar constraints f^o = 0.43^'^, 
with a slightly better reduced x 2 (x 2 — 115). The con- 
straint put on f^o for a flat Universe is remarkably consis- 
tent with the constrain derived from luminosity distances 
to SNI: Q Q = 0.3 ± 0.15 at the 90% confidence error only 
taking into account statistical errors (Fig. 7 in Perlmutter 
et al. 



19991) 



We check the robustness of our results on the self- 
similarity of clusters, with respect to the cluster temper- 
ature, by dividing the distant cluster sample in two equal 
sub-samples. We consider the favored cosmological model 
(Oo = 0.4, A = 0.6) and only data points with x > 0.08 
(corresponding to A^ p = 150 data points in total). We ob- 
tain a reduced x 2 of 0.61 (x 2 = 37 for 61 data points) 
for the kT < 6.6 keV subsample (13 clusters) and a re- 
duced x 2 of 0.88 (x 2 = 78 for 89 data points) for the 
kT > 6.6 keV subsample (12 clusters). Similarly, splitting 
the sample with respect to the cluster redshifts, yields a 
reduced x 2 of 0.82 (x 2 = 80 for 98 data points) for the 
z < 0.4 subsample (15 clusters) and a reduced x 2 01 0.67 
(x 2 = 35 for 52 data points) for the z > 0.4 subsample 
(10 clusters). In conclusion, for the favored cosmology, the 
distant cluster data for each individual subsample are in 
excellent agreement with the local reference profile. This 
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EM(x) oc A7 3 / 2 (l + z y 9/2 \s(9) {l + zf 
x cc 9 d A (z)A 1 J 2 (l + z) 3/2 



(20) 
(21) 



Scaling the observed and dimming corrected S(9) pro- 
file corresponds to translating it in a log-log plane. On 
the one hand, there is the translation of \og(dA(z)) along 
the x direction related to the conversion of angular radius 
into physical radius. On the other hand, the cluster cos- 
mological evolution requires an additional translation of 
-(3/2) log(A z ) - (9/2) log(l + z) in the y direction, and 
of (1/2) log(A z ) + (3/2) log(l + z) in the x direction, i.e 
along a line of slope -3. 

The cosmological parameters appear in the angular 
distance d\(z) and through the cluster over-density fac- 
tor A z . In the log- log plane, the scaled profiles for two dif- 
ferent cosmological models simply differ by translations. 
At a given redshift, varying the cosmological parameters 
simply corresponds to the same translation in the log-log 
plane for all the scaled profiles. The cosmological param- 
eters can thus only be constrained by comparing profiles 
at different redshifts. 

The reference scaled profile is determined from nearby 
cluster data. This profile depends itself on the cosmolog- 
ical parameters. At low redshifts, increasing fio (for both 
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reinforces the validity of the considered self-similar model, 
in particular the redshift dependence of the scaling rela- 
tions. 

An open model (A = 0) is also formally consistent 
with the data. However, the % 2 value keeps decreasing 
with decreasing fio, preventing a strict definition of the 
constraints. We thus only note that for fio = 0.1, we ob- 
tain a, x 2 — 78, similar to the value for the best model 
for the flat case. All open models with Slo > 0.17 give 
X 2 values which are larger than the values corresponding ^ 10 
to 90% confidence range of the flat case. Furthermore, we 
cannot consider arbitrarily low fio values. Obviously fio 
must be greater than the baryonic density derived from 
pri mord ial nucleosynthesis (fib = 0.03 — 0.06, Suzuki et 
al. E00l|) . In addition, the various approximations of the 
scaling models (in particular to compute the over-density) 
become less and less valid as fio decreases. 

5.3. Origin of the constraint on the cosmological 
parameters 

Comparing the scaled emission measure profiles of clus- 
ters at different redshifts appears to be a powerful method 
to constrain the cosmological parameters. To understand 
better the origin of the constraint, we examine in more 
details the variation of the scaled profiles, EM(x), with 
the cosmological parameters and redshift. 

It is useful to first explicitly identify this dependence, 
which is somewhat complex. The observed quantities are 
the surface brightness profiles S(0), which we correct for 
the (1 + z) 4 dimming factor. Combining Eq. [l^ and Eq. |l|, 
together with Eq. |^ and identifying the relevant factors, 
we can write: 
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Fig. 5. Dependence of the scaled profiles on the cosmolog- 
ical parameters assumed. Thick full line: reference scaled 
profile from the nearby cluster sample derived for fio = 1 . 
Dotted line: same for fio = 0.4, A = 0.6. At large radii, 
the logarithmic slope of the profile is ~ —3 (thin full line), 
corresponding to a /3-model with f3 = 2/3. Thick dashed 
line: Scaled profile of a z = 0.6 cluster one would derived 
assuming fio — 1 , if the 'true' cosmological model was 
fi = 0.4, A = 0.6, i.e if the 'true' scaled profile would fol- 
low the dotted line. Thin dashed line: same only taken into 
account intrinsic cosmological evolution. By using wrong 
values of the cosmological parameters, the derived profiles 
of distant and nearby clusters do not coincide any more. 
The discrepancy is essentially due to the dependence of 
the derived profiles on the angular distance assumed (see 
text for full discussion). 



flat and open models) mainly affects the A z factor, the an- 
gular distance being almost insensitive to the cosmological 
parameters. As A z increases with fio, increasing fio moves 
the scaled profile down and to the right, along the line of 
slope -3 (defining the scaling translation due to cosmolog- 
ical cluster evolution). This is illustrated in Fig. [|, where 
we compare the reference profile from the nearby cluster 
sample for two flat cosmological models (fio = 0.4, full 
line, and fio = 1, dotted line). A remarkable feature is 
the coincidence of the scaled profiles for the two models 
at large radii. This is due to the coincidence between the 
slope of the scaling translation (-3) and the slope of the 
profile at large radii (thin line in Fig. ||), so that the pro- 
file is just translated 'along it self. Note that this slope at 
large radii simply corresponds to a /3-model with (3 = 2/3, 
which was shown to fit well the mean profile of nearby 
clusters (Neumann & Arnaud 1999). At smaller radii, the 
slope of the profile becomes smaller than -3. As a result, 
the scaled profile for a high fio value always lies below the 
corresponding profile for a lower value of fio (Fig. ||). 

Let us now consider a high redshift cluster and assume 
that the correct cosmological model is a flat Universe with 
fi = 0.4. The scaled profile of this cluster, derived for 
fio = 0.4, will follow the corresponding reference curve. 
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However, this will not be the case if we assume another 
Qq value. A z varies more rapidly with Qq as the redshift 
increases. As a result, the scaled profile of a high redshift 
cluster is more affected by a change of fio than the scaled 
profile of a lower redshift cluster. Taking only into account 
the translation related to cosmological evolution, we com- 
pare the scaled profile of a z = 0.6 cluster, one would 
obtain assuming fio = 1 (thin dashed line in Fig ||) to the 
corresponding reference curve obtained for this cosmolog- 
ical parameter. The profiles still coincide at large radii, 
for the reason explained above. At small radii the z = 0.6 
profile is below the reference curve. It must be noted that 
this effect is small above 0.1 virial radius (less than 20%), 
i.e in the radial range considered. However, at high red- 
shifts, we have also to take into account the variation of 
the angular distance with Oo, which decreases with in- 
creasing Slo- The profile has to be further moved to the 
left, along the x axis, by \og(d A (z)^ 0=1) /d A (z) {no=OA) ^. 
At all radii, the profile of the z = 0.6 cluster (thick dashed 
line in Fig ||) does not coincide anymore with the reference 
profile. For a profile shape varying roughly as x~ 3 at large 
radii, the effect of d A (z) on the scaled emission measure 
at a given scaled radius is large, oc d A (z) 3 . At z = 0.6, the 
angular distance is about 18% higher for fio = 1 than for 
r^o = 0.4 and the profile is ~ 60% below the corresponding 
reference profile. 

This is exactly what we observe in Fig. ||. For Qq = 0.4 
(right panel) distant cluster profiles are consistent with 
the reference curve. For f2rj = 1 (left panel), all the data 
are moved down. The decrease is more important for dis- 
tant clusters, which are now systematically lower than the 
reference curve of nearby clusters, allowing us to exclude 
this cosmological model. The same reasoning applies for 
an open Universe. However, the variation of A z and d A (z) 
with Q,q is less pronounced for an open Universe than for a 
flat Universe. The differential effect is less important and 
the profiles coincide only for lower values of f2o- 

In conclusion, increasing f2rj moves the scaled profile 
of a given cluster down and right and decreases the scaled 
emission measure at a given scaled radius. However, in 
the radial range considered, the derived scaled profiles of 
any distant cluster, when compared to the corresponding 
reference profile, is mostly sensitive to the angular dis- 
tance d A (z) at the cluster redshift. This is this depen- 
dence, which essentially allow to constrain the cosmolog- 
ical parameters, via the well known dependence of d A (z) 
with fig and A. 



It is thus still below the profile corresponding to a lower 
f2o value. It can be also shown that in Eq. |2l] the product 
A^ 2 d A {z) (proportional to the inverse of the angular virial 
radius) always increases with increasing flo, and the profile 
remain globally moved to the right. The direction of the trans- 
lation, down and right, is readily apparent when comparin 
individual data points in the left and right top panels of Fig 



6. Evolution of the Lx-T relation 

We found that the scaled EM profiles of distant clusters 
coincide with the profile of nearby clusters, using a flat 
cosmological model with fio = 0.4. This means that the 
surface brightness profiles of distant clusters follow the 
evolution with redshift expected in the self-similar model, 
for this set of parameters. Since the X-ray luminosity is 
nothing else than the integral of the surface brightness 
profile, the evolution of the Lx~T relation should therefore 
also comply with this model. We check this point now. 

For consistency, we use the bolometric cluster luminos- 
ity estimated from the ROSAT data presented here, rather 
than ASCA. For each cosmological model, the total emis- 
sion measure within the virial radius is estimated by inte- 
grating the profiles up to the detection radius. The contri- 
bution beyond that radius was estimated using a /3-model 
with a slope /3 = 2/3, normalized to the emission measure 
at 0.3 virial radius. The luminosity was then estimated us- 
ing the cooling function computed with a MEKAL model 
at the cluster temperature (Tab. [l]). The ROSAT luminos- 
ity values, computed for Qq = 1 are given in Tab. |. They 
are in good agreement with the corresponding ASCA esti- 
mates from the literature. The median ratio between the 
two estimates is 0.97, with a standard deviation of 0.2 and 
there is no specific trend with redshift. 

The considered self-similar model assumes a non- 
evolving slope of the EM-T scaling relation, consistent 
with the slope (2.88) o f the l ocal Ly^-T relation established 



by Arnaud & Evrard ( 1999 ). We thus study the evolution 
of the normalisation of the L-&-T relation, assuming a con- 
stant slope of 2.88. For each cluster, we define, as in Sadat 
et al. (1998), the quantity: 



Cobs — 



Lo.o 5 T 12 ' 



(22) 



where Lx is the measured bolometric luminosity and Lo.05 
is the normalisation of the local Ly^-T relation, taken from 
the nearby sample (excluding A780, see below). This nor- 
malisation i s per fectly consistent with the data of Arnaud 
& Evrard ( 1999 ), the ratio of the two normalizations is 
0.99 (for Oq = !)• From Eq. [l7|, this quantity should evolve 
as: 



Cmod(-Z) 



A 



0.05 



1/2 



1 



1.05 



3/2 



(23) 



For consistency, we have normalized the theoretical func- 
tion, to the value at z = 0.05, the median redshift of 
the nearby sample. The observed C b s values, with error 
bars estimated from the errors on luminosity and temper- 
ature, are compared to the theoretical curve, Cmod(^)? in 
Fig. ^| for a critical density Universe (left panel) and for 
our best fit model (firj = 0.4, A = 0.6, right panel). Four 
clusters (A780, MS2137, MS1512 and RXJ1347) stand 
out with particularly high luminosities. Strong cooling 
flow clusters are known to lie above the Lx~T relation of 
weak or non-cooling flow clusters, considered by Arnaud 
& Evrard (|l999|). As mentioned above, MS2137, MS1512 
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Fig. 6. Evolution of the normalisation of the Lx~T relation for a critical universe (left panel) and a flat Universe 
with Qq — 0.4 (rtight panel). Data points: Observed bolometric luminosity divided by the luminosity estimated from 
the local Lx~T relation and the cluster temperature. Error bars include errors on the luminosity and temperature. 
Full line: theoretical evolution for the chosen cosmology and the self-similar model of cluster fo rmati on. Dotted line: 
the theoretical curve offset by the intrinsic scatter in log(Xx) estimated by Arnaud & Evrard ( 1999 ). The evolution 
of the normalisation of the Lx~T relation is consistent with the self-similar model for (fig = 0.4, A = 0.6). 



and RXJ1347 are indeed known to host strong cooling 
flows and the same is true for A780. They are thus dis- 
carded in the following. 

For a critical density Universe most of the data points 
lie below the theoretical curve, indicating that the distant 
clusters are under-luminous as compared to the theoreti- 
cal expectation. This is the same effect as observed for the 
scaled emission measure, which was also found to be too 
low as compared to the reference curve. We computed the 
X 2 value of the distant cluster data points about both the 
theoretical curve and a constant value of 1, correspond- 
ing to no evolution. We assume an in trinsic scatter of 0.13 
in iog(Lx) (Arnaud & Evrard 1999] ). We obtained a re- 
duced x 2 °f 2-3 and 1.1, respectively. We thus note that 
the observed C bs values seem actually consistent with no 
evolution at all, as found in previous studies of the Lx~T 
relation for fi = 1 ( Sadat et al. |l998j Fairley et al. |200C ; 
1999b . 



Reichart et al. 



On the other hand, the distant cluster data are consis- 
tent with the expected evolution of the Lx~T relation, for 
O = 0.4, A = 0.6. The x 2 , in this case, is \ 2 = 25 for 22 
clusters (reduced x 2 e( j = 1-1) • This good x 2 also indicates 
that the data are consistent with no evolution of the slope 
of the Lx~T relation with z, as we assume. By comparison, 
a reduced x 2 01 X 2 c d = 2-0 when the data are compared 
with the no evolution curve. The origin of the improve- 
ment, as compared to the flo = 1 model, is a combination 
of two factors i) the measured luminosity is higher for 
lower Oo, because the estimated distance is larger and ii) 
the expected evolution of the Lx~T relation is more mod- 
est due the factor A z , a decreasing function of z. We find 
that C mo d(z), for this cosmological model, is well approx- 
imated by a power law over the redshift range considered 
here (z = 0. - 0.83): C mod (z) = (1 + z)^, with r) = 1.23. 



This can be compared with the results of Reichart et al. 
flLggg ), who assume such a dependence. Their results are 
based on a compilation of data from the literature for clus- 
ters with z < 0.5. For fl = 0.4, A = 0.6 (q = -0.4) they 
found r) = 1.36^° their best fit is thus also in good 
agreement with the theoretical expectation. 

In summary, as expected from our study of the emis- 
sion measure profiles, the normalisation of the Lx~T re- 
lation, derived for a flat low density Universe, does evolve 
with redshift. The observed evolution is consistent with 
the self-similar model of cluster formation. 



7. Discussion 

7.1. Cluster formation and evolution 

We found an excellent agreement, both in shape and nor- 
malisation, of the set of scaled profiles of distant clusters 
with the reference nearby scaled, for a low density flat 
Universe. This indicates that hot galaxy clusters consti- 
tute a homologous family up to high redshift, with the 
cluster properties scaling with z as expected in the sim- 
plest self-similar model. These scaling laws are derived 
from the assumption that clusters form at fixed den- 
sity contrast as compared to the critical density of the 
Universe. Our results thus support this standard picture 
for the gravitational collapse of the dark matter compo- 
nent. 

Consistently, the evolution of the normalisation of the 
Lx~T relation was found to comply with the self-similar 
model for O = 0.4, A = 0.6. The apparent inconsistency 
with this model, claimed in some previous studies, was in 
fact due to the a priori choice of a particular cosmological 
model (Qq = 1), and is not, per se, an indication of extra- 
physics. That does not mean that such physics does not 
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exist. We emphasize again that the present study only con- 
cerns relatively hot clusters, for which non-gravitational 
effects like pre-heating is minimal, and that the slopes of 
the M gas -T and Lyj-T relations remain inconsistent with 
the standard scaling laws. The simplest self-similar model 
is clearly insufficient to fully describe the properties of 
the gas component of the ICM. Our data are consistent 
with no evolution of these slopes with z but much better 
temperature estimates and larger samples are required for 
further tests of this point. If this was confirmed, the em- 
pirical slopes, together with similarity in shape up to high 
z, have to be explained in terms of the specific gas physics 
in the picture of structure formation. 



7.2. A new method to constrain the cosmological 
parameters from clusters? 

The basic assumption of the method to constrain the cos- 
mological parameters, that we validate a-posteriori, is that 
clusters form a homologous population. In particular, once 
scaled according to cluster temperature and redshift, the 
EM profiles, derived from the observed surface brightness 
profiles S(9), follow a universal function, common to all 
clusters. This universal function, determined from nearby 
cluster data, is then used as 'standard candle' to constrain 
the cosmological parameters. The scaled profiles are not 
used as simple distance indicators, in the strict sense. They 
do depend on the angular distance, via the conversion of 
angular radius into physical radius, but, in the self-similar 
model, the evolution of the scaling relations used in the 
scaling process also depends on the cosmological param- 
eters. However, we have shown that the second effect is 
actually negligible. The scaled profiles show a strong de- 
pendence on the angular distance, oc <Ia(z) 3 at any given 
scaled radius. This makes them competitive with other 
distance indicators like SNI (ex cIa(z) 2 ) and indeed, we 
were able to derive constraints of similar quality. 

The method is however by no mean as powerful as 
the cosmological test based on the evolution of the clus- 
ter mass function, N(M,z), which depends much more 
strongly on f2o- For instance, the abundance of massive 
clusters falls by an order of magnitude by z ~ 0.6 in a 
critical density Universe, while it remains unchanged in an 
fto = 0.3 Universe (Blanchard & Bartlett 1998 ). However, 
indirect mass indicators, like T or L%, have to be used 
and the selection function for flux-limited cluster surveys 
depends on the cluster scaling and structural properties. 
This N(M, z) test, by nature, thus requires a good knowl- 
edge of the cluster scaling properties, and their evolution 
with z, that we study here. The major advantage of the 
present method is that it is more direct, i.e directly based 
on observed quantities. 

The proposed method suffers, nevertheless, from in- 
trinsic limitations. First, there is some intrinsic dispersion 
in the cluster properties. The typical dispersion is about 
±25%, which is not so different from the ~ 60% in cIa{z) 3 
between an fio = 0.4 (A = 0.6) and an Qq = 1 cosmology 



for a cluster at a redshift as high as z = 0.6. Therefore, 
the method requires to consider a large sample of clusters, 
the scaled profiles of distant clusters will coincide with the 
reference profile, for the correct values of the cosmologi- 
cal parameters, only in a statistical sense (i.e. on average) 
and not on a cluster by cluster basis. However, this intrin- 
sic dispersion can be measured and explicitly taken into 
account. This is done in the present analysis, where the 
set of distant cluster profiles is compared to the reference 
profile, with its dispersion, using \ 2 statistics. 



A more serious concern is that the method is intrin- 
sically model dependent. As for other distance indica- 
tors, we cannot exclude that some evolutionary effects, 
not taken into account, bias the results. In particular, the 
present analysis relies on the assumption that the slope of 
both the EM-T and Ry-T relations does not evolve with 
z. If this was not the case, different constrains on might 
be obtained. For instance, a decrease of the EM-T rela- 
tion slope with z would boost the scaled emission measure 
profiles of high z clusters, as compared to low z, and could 
possibly mimic the effect of a low f2o Universe. To further 
quantify this point, we let the slope of the EM-T relation 
vary with redshift. As a test case, we consider a very sim- 
ple model. With respect to the standard EM-T relation, 
EM cx Az /2 (1 + zf^T 1 - 38 , the EM-T relations at the 
various redshifts are rotated, around a common reference 
point of given temperature T p , so that the slope varies 
linearly with z: EM' = EM(T/T p )-( 5a< - z » . The slope at 
z = 0.05, the median redshift of the nearby cluster sam- 
ple is kept unchanged and we consider a slope change as 
high as ±20% at z = 0.8: 8a{z) = ±0.37(z - 0.05), corre- 
sponding to EM cx T 11 or EM cx T 166 at z = 0.8. The 
scaled emission measure profiles are then given by Eq. [l9[ 
modified by a multiplicative factor of (T /T p )~ Sa( - z \ The 
change in the EM-T relation, and thus its impact on 
the fio estimate, depends somewhat critically on T p ; it 
is obviously maximal for low values of T p . If we take as 
reference temperature, T p = 3.5 keV, the minimum tem- 
perature of the sample, the scaled emission measure of a 
10 keV cluster at z — 0.8 is increased by about 35%, for 
a 20% decrease in slope. As expected, higher f2 values 
are then derived: Q — ®-6-o.i5 ( nat Universe). If instead 
we assume a steepening of the EM-T relation with z, 
we obtain f^o = 0.25 ± 0.11. Both evolution models are 
consistent with the data, the corresponding reduced \ 2 is 
not significantly different from the one obtained for the 
standard model. We note however that, in both cases, a 
critical density Universe remains excluded. Furthermore, 



when we take T p = 6.5 keV, the median temperature of 
the sample, the slope evolution has a negligible effect on 
the derived parameters: we obtain Slo = 0.47 ± 0.14 and 
f2 = 0.34+q ^ for a +/ — 20% change in slope, respec- 
tively. This systematic effect is twice smaller than the sta- 
tistical errors on the f^o value. Finally, we also checked the 
effect of a similar change in the slope of the Ry-T relation. 
The effect is similarly negligible: we obtain f^o = 0.30lg 09 
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and Oo = 0.471q 13 f° r a 20% increase and decrease of the 
slope respectively. 

Our conclusions also depend on the assumed EM-T 
relation, because we are not considering clusters of similar 
temperatures at all redshifts. The median temperature is 
5.8 keV for the nearby sample, and 6.6 keV and 7.9 keV for 
distant clusters in the redshift ranges 0.3-0.43 and 0.45- 
0.83, respectively. For a shallower EM-T relation, like the 
standard scaling relation, most of the scaled profiles of dis- 
tant clusters, lie above the scaled profiles of nearby cluster 
for Q = 0.3 (Fig. |). In that better agreement be- 

tween the distant cluster data and the nearby cluster data 
is obtained for 0,$ = 1 than for Slo = 0.3. However, we 
emphasize again that the standard scaling law is not con- 
sistent with the slope of the Lx—T relation and that the 
adopted scaling relation decreases significantly the scatter 
of the scaled profiles. 

Nevertheless, a precise determination of the scaling 
with temperature, and of its possible evolution, remains 
essential to achieve a fully consistent description of cluster 
evolution and to assess possible systematic errors on the 
cosmological parameters. Again, the quality of the present 
data is rather insufficient for high precision tests. Finally, 
there will always remain the possibility of a degeneracy 
between the evolution of the cluster properties and the 
variation of the angular distance with redshift. Such sys- 
tematic errors can only be assessed by comparing the re- 
sults obtained by various methods. The good agreement 
obtained between our results and the constraints based 
on the luminosity distance of SNI is an encouraging sign 
that both methods are unbiased and the underlying mod- 
els correct. 

7.3. Comparison with previous work based on the 
size-temperature relation 

The cosmological parameters can also be constrained, as 
proposed by Mohr et al. (200C), using the isophotal size - 
temperature (ST) relation. Mohr et al. (2000) showed that 
the normalisation of this relation is insensitive to cluster 
cosmological evolution, considering the same model than 
in the present study. Their test of the cosmological pa- 
rameters is thus made via the dependence of the size on 
the angular distance. Both this method and ours thus use 
quantities derived from cluster surface brightness profiles, 
as distance indicators. 

The main difference between the two methods is that 
we consider scaled quantities rather than physical quanti- 



ties. At large radii, considered by Mohr et al. (2000), this 
is equivalent. Due to the coincidence between the slope 
of the scaling translation related to cluster cosmological 
evolution and the slope of the profiles (see Sect. 5.3), the 
profiles of all clusters coincide at large radii, both in the 
scaled space and in the physical space. This is the origin 
of the invariance of the isophotal size with redshift (the 
arguments developed by Mohr et al. ( [2000 ) are actually 
similar). For an isophotal size evaluated from cluster im- 



ages, the method of Mohr et al. (200C) is equivalent, in our 
approach, to consider only data points at a given scaled 
emission measure. 

Our method, where we consider the whole set of data 
points, can thus be regarded as a generalization of the 
method proposed by these authors. Note that, by work- 
ing with scaled quantities, we are able to consider data 
at small radii, where cosmological evolution has to be 
taken into account (even if it does not depend sensitively 
on the cosmological parameters). The method we pro- 
pose presents several advantages. Obviously tighter con- 
strains can be obtained by considering the whole set of 
data points. No parametric fit of the surface b rightn ess 
profiles, as the one introduced by Mohr et al. ( 2000 ), is 
required. Furthermore, it allows a more complete test of 
the underlying self-similar model. 

We stress on the agreement between the results ob- 
tained Mohr et al. ( 2000| ) and ours. Both studies exclude 
a critical density Universe. Low Qq values are favored, 
somewhat lower when we use the evolution of the ST re- 
lation. One also notes that lower Oo values are favored for 
an open model than for a flat model. 

8. Conclusion 

In this work based on ROSAT data and published ASCA 
temperatures we study the surface brightness profiles of a 
sample of hot (kT > 3.5 keV) galaxy clusters, covering a 
redshift range z = 0.04 — 0.83. For both open and flat cos- 
mological models, the derived emission measure profiles 
are scaled according to the self-similar model of cluster 
formation. We use the standard scaling relations of cluster 
properties with redshift. The physical radius is normalized 
to the virial radius, estimated from the classical virial re- 
lation. The slope of the EM-T relation depends on the 
assumed slope of the M gas -T relation. We consider both 
the standard scaling relation M gas cx T 15 and the em- 
pirical local relation M gas cx T 1 - 94 (Neumann & Arnaud 



( [200 1| ), assuming the slope does not evolve with z. 

Our analysis of the scatter of the scaled profiles, sug- 
gests that the empirical slope of the M gas -T relation fits 
better the cluster properties than the standard value, over 
the whole redshift range z — 0.04—0.8. As for nearby 
clusters, a large dispersion in the central core is observed, 
and we therefore consider only the region above typically 
0.1i? v - 

Applying the empirical EM-T relation, the set of 
scaled profiles of the distant cluster sample are compared 
to the average scaled profile of nearby clusters, using a \ 2 
test. An excellent agreement, both in shape and normalisa- 
tion, of the distant cluster data with this reference nearby 
scaled profile is obtained for a flat low density Universe 
(see also below). Consistently, the evolution of the normal- 
isation of the Lx~T relation was found to comply with the 
self-similar model. The apparent inconsistency with this 
model, claimed in some previous studies, was in fact due 
to the a priori choice of a particular cosmological model 

(n„ = 1). 
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This indicates that hot galaxy clusters constitute a 
homologous family up to high redshifts and supports 
the standard picture for the gravitational collapse of the 
dark matter component. However, the simplest self-similar 
model is insufficient to fully describe the properties of 
the gas component of the ICM, as indicated by the non- 
standard slope of the M gas -T (and Lx~T ) relation. If 
confirmed, this slope, together with similarity in shape up 
to high z, have to be explained in terms of the specific gas 
physics in structure formation scenario. 

Because of the intrinsic regularity of the hot cluster 
population, we showed that the scaled emission measure 
profile, determined from nearby cluster data, can be used 
as 'standard candle' to constrain the cosmological param- 
eters, the correct cosmology being the one for which the 
profiles at different redshifts coincide. The scaled profiles 
of distant clusters, as compared to the reference profile, 
mostly depend on the angular distance, as g?a(.z) 3 , mak- 
ing them powerful distance indicators. The method is, in 
addition, more powerful than the test based on the size- 
temperature relation (Mohr et al. 200C ), because it uti- 
lizes the full information contained in the cluster profiles, 
rather than a particular point of the profiles. 

Using this new method, we were able to exclude a 
critical-density model (f2 = 1) (at 98% confidence level). 
The data favor a flat Universe with a low matter density, 
even if the open model is not formally excluded. We find 
a value of fl Q = 0.40±g;^ (at 90% confidence level). This 
test relies on the fact that we are using the right scal- 
ing relations, in particular for the EM-T relation. It is 
thus, by nature, a model dependent method, although the 
model can be to some extent, validated a posteriori. 

At this stage, our proposed method has to be taken 
more in terms of an independent consistency check of 
the constraints on cosmological parameters rather than 
"an ultimate cosmological test". The constraint derived 
on f2o is in remarkable agreement with the constraint ob- 
tained from luminosity distances to SNI (Perlmutter et al. 



199E) or from combined analysis of the power spectrum of 



the 2dF galaxy redshift Survey and the CMB anisotropy 
(Efstathiou et al |2002| ; see also Melchiorri et al |200C ; 
Stompor et al. 200C; de Bernardis et al. 2002; Pryke 



et al. 2001 and references therein). This is an additional 
sign that we are entering an era where cosmological tests 
converge and we can expect that soon the cosmological pa- 
rameters will be accurately known. In this context, cluster 
scaling and structural properties will be more adapted to 
test the physical processes in the structure formation pic- 
ture. Significant progresses in this field require high qual- 
ity data with measurements down to the virial radius that 
will be provided by the new generation of X-ray observa- 
tories (Chandra and XMM-Newton) . They also require a 
large sample of distant and nearby clusters so that i) the 
intrinsic dispersion is pinned down, ii) we improve our 
knowledge of the local relations and the temperature and 
dark matter profiles, and iii) we fully assess the evolution 
with z. 
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Appendix A: Analytical fit of the reference scaled 
profiles 

For a critical Universe the reference profile derived from the 
set of nearby cluster profiles is well fitted by a /3-model above 
a scaled radius of x = 0.1: 



EM(x) = EMq 



1+1^ 



-3/5+0.5 



(A.l) 



with EM = 0.945, x c = 0.116 and (3 = 0.67. This analytical 
formula is accurate to ~ 7% (with typical errors less than 5%) 
in the range x = 0.1-0.7. We emphasize that this formula must 
not be used for lower values of x (where the scatter increases 
and the profiles are on average more peaked), as well as above 
x = 0.7, corresponding to the last measured point. 

The individual cluster scaled profiles, derived from the ob- 
served surface brightness profiles, depend on the cosmologi- 
cal parameters, via the factor A z and the angular distance, 
as given in Eq. ^l|. The reference profile, for any cosmologi- 
cal model, is well approximated by simply scaling the Qg — 1 
reference profile, with the A z and d\{z) factors estimated at 
z — 0.05, the mean redshift of th e sam ple. The reference profile 
is thus given by a /3-model (Eq. Al) with: 



EM = 0.945 A~ 



3/2 



,V2 



0.67 



(A.2) 
(A.3) 
(A.4) 



where /d A is the angular distance normalised to its value for 
fio = I- The A0.05 factor is given Eq. | and Eq. g with z — 
0.05. It is accurately (within less than 0.2%) aproximated by 
the following polynomial expression: 

Ao.os = 1 - 0.365(^0 - 1) - 0.174(fi - l) 2 for fi < 1, A = 
Ao.os = 1 - 0.522(ft () - 1) - 0.199(fio - I) 2 for fi + A = 1 

(A.5) 

Similarly the /d A factor can be aproximated by: 

,fd A = 1 + 1.25 10~ 2 (ft - 1) for n < 1, A = 

/d A = 1 + 3.75 10~ 2 (fi - 1) for fl + A = 1 (A.6) 

The overall accuracy is similar to the accuracy obtained 
for the critical Universe reference profile, for Qo > 0.1. Again 
we emphasize that the analytical formula must only be used 
between x min = 0.1Aj!/ 2 5 /d A and x max = 0.7Aj' 2 6 / dA . 
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Appendix B: Computation of the x 2 value 

Here, we describe the way we compute the x value used to 
compare the set of scaled emission measure profiles of distant 
clusters to the reference profile. 

The data consist of a set of scaled emission measure, yk = 
EMk, measured at the scaled radius Xk- These quantities are 
derived from the surface brightness, SBk, at angular radius 9k 
(corresponding to the scaled radius Xk) and the temperature 
of the specific cluster, Ty. 



y k oc SB k T k 



Xk 



oc T, 



1/2 



(B.l) 
(B.2) 



where a is the slope of the EM-T relation. The corresponding 
errors, a Vk and a Xk are: 



vl 



1 1 

SBl +a Tl 



(B.3) 
(B.4) 



where <jT k and (?sB k are the uncertainties on Tk and SBk re- 
spectively. The errors on Xk and yk are thus correlated through 
the error on the temperature. The correlation factor pk is: 



PkCx 



x k yk 



a 
2 



(B.5) 



Let us note Y = f(X) the equation of the reference curve 
to which this data set is compared. In practice, it is given in 
tabular form and the data for any value of X is obtain ed by 
interpolation. The \ 2 expression can be found in York (L969) 
for the case of correlated errors: 



(B.6) 



where N is the number of data points and Sk the distance of 
the data point k to the reference curve, which is obtained by 
minimizing over X the function: 



Sk(X) 



1-Pl 



(xk-X) 2 , (y k -f(X)) 2 



-2p, 



+ 



(x k - X) (y k - f(X)) 



(B.7) 



Since the reference function is not linear, this minimization, 
which actually determines the 'closest' reference point, is done 
numerically. Up to this stage, we have not taken into account 
the dispersion, (Tftx) a t radius X, observed around the refer- 
ence curve. This is done by a dding quadratically this dispersion 
to a Vk . Eq. |B,5| and Eq. B.7 remain the same, where a Vk is re 

placed by \j(J% 
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